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Lorentz violation at high energies might lead to non linear dispersion relations for 
the fundamental particles. We analyze observational constraints on these without 
assuming any a priori equality between the coefficients determining the amount 
of Lorentz violation for different particle species. We focus on constraints from 
three high energy processes involving photons and electrons: photon decay, photo- 
production of electron-positron pairs, and vacuum Cerenkov radiation. We find 
that cubic momentum terms in the dispersion relations are strongly constrained. 
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1 Introduction 

There are several reasons to suspect that Lorentz invariance may be only a 
low energy symmetry. This possibility is suggested by the divergences of local 
quantum field theory, as well as by tentative results in various approaches to 
quantum gravity. Moreover, the unboundedness of the boost parameter makes 
experimental verification of exact Lorentz symmetry impossible in principle. 

One can study the possibility of Lorentz violation, without a particular 
fundamental theory in hand, by considering its manifestation in dispersion 
relations for matter. It is natural to assume that such dispersion relations 
E 2 (p) can be characterized by an expansion with integral powers of momentum, 

E 2 = p 2 + m 2 + [Ap 2 + Bp 3 /K + Cp 4 /K 2 + 0(p 5 )] . (1) 

Here A, B, C are dimensionless coefficients which might be positive as well 
as negative and Kq is the "quantum gravity" scale (often identified with the 
inverse Planck length, Kq = 2ir/Lp{). [Throughout this paper p denotes the 
absolute value of the 3-momentum vector p, and we use units with the low 
energy speed of light in vacuum equal to unity] 

Different approaches to quantum gravity suggest different leading order 
Lorentz violating terms. The terms with coefficients A, B, C have mostly been 
considered so far. Since the p 2 term is not suppressed by the Planck scale, it 
might be thought to be largest, however observations severely limit A to be 
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much less than unity (see e.g. a and references therein). The higher order terms 
are naturally small, with coefficients of order unity. 

Our strategy here is to take a purely phenomcnological stance and consider 
the constraints that high energy observations impose on dispersion relations of 
the form 

El^pl+ml + riapl/Kr 2 , (2) 

where a labels different fields and n > 3. In the absence of a fundamental the- 
ory one has no reason to expect any particular relation between the coefficients 
r\ a for different particles, except perhaps that they should all be of the same 
order of magnitude. 

Dispersion relations of the kind (|^) produce kinematic constraints from 
energy-momentum conservation that differ from the usual Lorentz invariant 
case. As a result reactions can take place that are normally forbidden, and 
thresholds for reactions are modified. Observational consequences may seem 
out of reach because of the Planck scale suppression of the Lorentz violating 
terms (assuming that, as generally expected, Kq is of order the Planck scale). 
However this is not so. One can expect deviations from standard kinematics 
when the the last two terms of (||) are of comparable magnitude.F] Assuming 
77 is of order unity this yields the condition pdev ~ ( m2 Ko~ Y > which is 
- (m/m e ) 2/3 x 10 TeV for n = 3 and ~ (m/m e ) 1 '' 2 x 10 4 TeV for n = 4. 
Although these energies are currently not achievable in particle accelerators 
(except in the case of the massive neutrinos which however are weakly coupled) 
they are in the range of current astrophysical observatdnns 

In fact, it has been suggested by several authors H'BHO'BBIj (see alsoB and 
references therein) that we may already be observing deviations from Lorentz 
invariance via the existence of two puzzles in modern astrophysics: the missing 
GZK cut-off on cosmic ray protons with ultra high energy greater than 7 x 10 19 
eV, and the apparent-overabundance of gamma rays above 10 TeV from the BL 
Lac system Mkr 501 Here we shall mostly not consider the constraints im- 
posed by asking Lorentz violation to explain these puzzles. Instead we restrict 
our attention to constrainta-jimposed by consistency with known phenomena 
(or lack thereof). (See alsoUBu for a similar discussion.) 

2 Observational constraints 

To find the strongest observational constraints without assuming a priori re- 
lations between the coefficients ij a of (||) we focus here on processes involving 

a In general of course one must look at the specific reaction in order to estimate the energies 
at which deviations from standard behavior can be expected. For example in the case of the 
reaction 77 — > e+e — it will be the electron mass that sets the scale. 
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just two fundamental particles, photons and electrons. We also restrict to the 
case n = 3, since this should be most tightly constrained. If it can be ruled 
out, one can then move on to the n = 4 case. 

The modified dispersion relations for photons and electrons in general allow 
two normally forbidden interactions: photon decay, 7 — » e + e~, and vacuum 
Cerenkov radiation, e~ — > e~ "7. If allowed these processes happen very rapidly. 
In addition the threshold for photon annihilation, 77 — > e + e~, is shifted. We 
consider constraints that follow from three obserjzations: (i) Gamma rays up 
to ~ 50 TeV of cosmological origin arrive on earthEj, so photon decay does not 
occur up to this energy, (ii) Electrons of energy ~ 100 TeV are believed to pro- 
duce observed X-ray synchrotron radiation coming from supernova remnants. 
Assuming these electrons are actually present, vacuum Cerenkov radiation 
must not occur up to that energy, (iii) Cosmic gamma rays below 10 TeV are 
believed to be absorbed in a manner consistent with photon annihilation off 
the IR background with the standard threshold. Observation (iii) is not model 
independent, so the corresponding constraint is tentative and subject to future 
verification. 

To derive the observational constraints one needs to determine the thresh- 
old for each the process, i.e. the lowest energy for which the process occurs. 
The details concerning our determination of the thresholds are reported else- 
whereliil. Assuming monotonicity of E(p) (for the relevant momenta p <C Kq) 
we have shown that all thresholds for processes with two particle final states 
occur when the final momenta are parallel. Moreover for two particle initial 
states the incoming momenta are antiparallel. This geometry has been assumed 
in previous works but to our knowledge it was not shown to be necessary. (In 
fact it is not necessary if E(p) is not monotonic.) 

To eliminate the subscript a we introduce £ ■= V-y an d V '•= T) e . The 
constraints will restrict the allowed region of the plane. For the rest of the 
paper we assume Kq is the Planck energy and we use units with Kq = 1. 

2.1 Photon decay 

Photon decay is allowed only above a broken line in the 77— £ plane given by 
£ = 7]/ 2 in the quadrant £, 77 > and by £ = 77 in the quadrant £, r\ < 0. Above 
this line, the threshold is given by 




for 77 < £ < 0. 



for £ > 0, 



(3) 



(4) 
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The first relation (|^) arises when the electron and positron momenta are equal 
at threshold. In standard Lorentz invariant kinematics such "equipartition" of 
momentum always holds at threshold for pair production. In all previous work 
on Lorentz violating dispersion it has been assumed to hold. Surprisingly, 
however, in the present case the threshold may occur with an asymmetric 
distribution of momentum. The second relation ^ applies in those cases. 
The constraint we impose is that the threshold is above 50 TeV, the energy of 
highest observed gamma rays from the Crab nebula E3. 

2.2 Photon annihilation 

The threshold relations for a gamma ray to annihilate with an IR background 
photon of energy e take approximately the same form as for photon decay, with 
the replacement ^ — > ^ , where £ = £ + Ae/k 2 h . (Here we have used the fact 
that e is much smaller than any other scale in the problem.) The two different 
relations arise for the same reason as in the case of photon decay, and they 
correspond respectively to cubic and quartic polynomial equations for fc t h- 

For the observational consequences it is important to recognize that the 
threshold shifts are much more significant at higher energies than at lower 
energies. To exhibit this dependence, it is simplest to fix a gamma ray energy 
k and to solve for the corresponding soft photon threshold energy e t h- Taking 
the ratio with the usual threshold eth.Oj we nn d a dependence on k at least as 
strong as k 3 / 2 . Introducing k w := fc/(10TeV), we have 

1 + 0.05 (r) - 2£) kf for £' > 0, (5) 

0.1 fa -0 k 3 w + -0.2 nk 3 10 for ry < £' < 0. (6) 

High energy TeV gamma rays from the blazars Markarian 421 and Markar- 
ian 501 have been detected out to 17 TeV and 24 TeV respectively. Although 
the sources are not well understood, and the intergalactic IR background is 
also not fully known, detailed modeling shows that the data are consistent with 
some absorption by photon annihilation off the IR background (see e.g.ll3'E3EJ 
and references therein). However, while the inferred source spectrum for Mrk 
501 is consistent with expectations for energies less than around 10 TeV, above 
this energy there are_£ar more photons than expected according to some IR 
background models Bta. 

The uncertainty in the blazar and IR background models currently pre- 
cludes sharp constraints from photon annihilation. Instead, we just determine 
the range of parameters £,77 for which the threshold fc t h lies between 10 TeV 
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and 20 TeV for an IR photon with which a 10 TeV photon would normally be 
at threshold. We choose this range since (i) lowering the threshold would make 
the overabundance problem worse, and (ii) raising this threshold by a factor 
of two might explain the potential overabundance of photons over 10 TeVa. 
If the overabundance puzzle exists and is to be resolved in this way, the effect 
of the threshold shifts must be enhanced for photons above 10 TeV relative 
to those below 10 TeV (since the latter seem to be well accounted for). This 
enhancement could arise partly from the shape of the IR background spec- 
trum, but the shift itself is also enhanced, as seen in the energy dependence of 
equations . The hypothesis BM that the potential overabundance is due 
to Lorentz violating dispersion with n — 3 therefore appears to be consistent 
with current observations. 



2.3 Vacuum Cerenkov radiation 

An electron can emit Cerenkov radiation in the vacuum if r\ > or if r\ < 
and £ < r). In this case there are also two threshold relations, depending on 
whether the threshold occurs with emission of a zero-energy photon or with 
emission of a finite energy photon. These two cases correspond to the two 
following relations, respectively: 

/ 2 \ V3 



Pth = [j^J for V > and £ - -3? ?' ( 7 ) 

1 /3 

(4 m? (£ 4- r? j \ 
^ for £ < -3?7 < or £ < 77 < 0. (8) 
(£ - V) J 

The vacuum Cerenkov process is extremely efficient, leading to an energy loss 
rate that goes like E 2 well above threshold. Thus any electron known to prop- 
agate must lie below the threshold. Electrons of energy ~ 100 TeV are believed 
to produce observed X-ray synchrotron radiation coming from supernova rem- 
nants. Thus for example in the region of the parameter plane where (|7|) holds 
we obtain the constraint r\ < m 2 /2p^ h ~ 10 -3 . 



3 Combined constraints 

Putting together all the constraints and potential constraints we obtain the 
allowed region in the 77-£ plane (see Figure l]) . The photon decay and Cerenkov 
constraints exclude the horizontally and vertically shaded regions, respectively. 
The allowed region lies in the lower left quadrant, except for an exceedingly 
small sliver near the origin with < r) < 10~ 3 and a small triangular region 
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(-0.16 < rj < 0, < S, < 0.08) in the upper left quadrant. The range of 



the photon annihilation threshold discussed in subsection 2.2 falls between 
the two roughly parallel diagonal lines. This intersects the otherwise allowed 
region in a finite, narrow wedge where £ and r\ are negative and of order unity 
(apart from a minuscule invisible region near the origin with r\ > 0). If future 




Figure 1: Combined constraints on the photon and electron parameters for the case n = 3. 
The regions excluded by the photon decay and Cerenkov constraints are lined horizontally in 
blue and vertically in red respectively. The photon annihilation threshold limits fall between 
the two diagonal green lines. The dashed line is £ = r\. 



observations of the blazar fluxes and the IR background confirm agreement 
with standard Lorentz invariant kinematics, the region allowed by the photon 
annihilation constraint will be squeezed toward the upper line (fc t n ~ k s ). If the 
overabundance of gamma rays from Mrk 501 is indeed due to Lorentz violation 
of the sort we are considering, the region will be squeezed toward the lower 
line (k t h ~ 2fc s ) or may even shift toward lower values of £ and 77 if a yet larger 
threshold shift is indicated. 
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4 Conclusions 



We have seen that a conservative interpretation of observations puts strong 
constraints on the possibility of Planck scale cubic modifications to the electron 
and photon dispersion relations. The alkiwed region includes £ = 77 = — 1, 
which has been a focus of previous workB'tm. The negative quadrant has most 
of the allowed parameter range. It is interesting to note that in this quadrant 
all group velocities are less than the low energy speed of light. 

To completely rule out the cubic case would require new observations. 
Finding higher energy electrons would not help much, while finding higher en- 
ergy undecayed photons would squeeze the allowed region onto the line £ = r\. 
To further shrink the allowed segment of this line would require observations 
allowing the usual threshold for photon annihilation to be confirmed to higher 
precision. Perhaps other processes could be used as well. If a priori relations 
among the coefficients in the dispersion relations for different particles are hy- 
pothesized, stronger constraints can of course be obtained. 

This research was supported in part by NSF grant PHY-9800967. 



References 



S. Coleman and S. L. Glashow, Phys. Rev. D 59, 116008 (1999). 
G. Amelino-Camelia and T. Piran, Phys. Rev. D64, 036005 (2001). 
L. Gonzalez-Mestres, "Gamma and cosmic ray tests of special relativity' 1 



astro-ph/0011181 and references therein. 



4. O. Bertolami, "Ultra-high energy cosmic rays and symmetries of space- 



time", astro-ph/0012462 



5. T. Kifune, Astrophys. J. 518, L21 (1999). 

6. W. Kluzniak, Astropart. Phys. 11, 117 (1999). 

7. R. Aloisio, et ai, Phys. Rev. D62, 053010 (2000). 

8. G. Sigl, Led. Notes Phys. 556, 259 (2000). 

9. R. J. Protheroe and H. Meyer, Phys. Lett. B493, 1 (2000). 

10. T. Tanimori et ai, Ap J 492, L33 (1998). 

11. T. Jacobson, S. Liberati, and D. Mattingly, Ln Preparation. 

12. F. Krennrich et al., "Cutoff in the TeV energy spectrum of Markarian 
421 during strong flares in 2001 ", astro-ph/0107115 . 

13. O. C. de Jager and F. W. Stecker, "Extragalactic gamma-ray absorption 
and the intrinsic spectrum of Mkn 501 during the 1977 flare," 



ph/0107103 



astro- 



14. F. A. Aharonian, A. N. Timokhin and A. V. Plyash eshnikov, "On the 
origin of highest energy gamma-rays from Mkn-501," astro-ph/0108419 . 



7 



